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ABSTRACT 

The contribution of higher-order multipoles to radiative and non-radiative decay of a 
single dipole emitter close to a spherical metallic nanoparticle is re-examined. Taking a 
Ag spherical nanoparticle (AgNP) with the radius of 5 nm as an example, a significant 
contribution (between 50% and 101% of the total value) of higher-order multipoles to 
non-radiative rates is found even at the emitter distance of 5 nm from the AgNP sur- 
face. On the other hand, the higher-order multipole contribution to radiative rates is 
negligible. Consequently, a dipole-dipole approximation can yield only an upper bound 
on the apparent quantum yield. In contrast, the non-radiative rates calculated with the 
quasistatic Gersten and Nitzan method are found to be in much better agreement with 
exact electrodynamic results. Finally, the size corrected metal dielectric function is shown 
to decrease the non-radiative rates near the dipolar surface plasmon resonance. 

PACS numbers: 33.50.-j, 33.50.Hv, 32.50. +d, 78.67.-n, 71.45.Gm, 73.22.Lp, 78.67.Bf 
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1 Introduction 



Metallic surfaces are known to increase the local electromagnetic fields and modify both 
excitation and emission rates of proximate fluorophores, chromophores, and quantum dots 
[USElllllSlElCOIHllHl ^3 ^0 02] ■ Especially promising is modifying fluorescence rates 
by means of metallic nanoparticles (MNPs) for addressing various issues of interest to 
biology [21 |3l HI [5]. The MNPs can quench fluorescence as much as 100 times better 
than other quenchers of fluorescence, such as DABCYL, and open new perspectives in 
the use of hybrid materials as sensitive probes in fluorescence-based detection assays 
[2]. Molecular beacons comprising the MNP acceptors can detect minute amounts of 
oligonucleotide target sequences in a pool of random sequences and provide an improved 
detection of a single mismatch in a competitive hybridization assay for DNA mismatch 
detection [2]. The MNP acceptors can also be employed for probing changes in distances 
for protein interactions on DNA using a molecular ruler approach called nanometal surface 
energy transfer [9]. The MNPs induced decrease in lifetime of fluorophores can result in 
an effective increase in photostability. All that carries potential for the next generation 
superstructures in clinical diagnostic, DNA sequencing, genomics, and biological detection. 
The latter is the subject of a new field of radiative decay engineering that aims at a 
predetermined modification of decay rates [31 E]. 

In theory, fluorophores, chromophores, and quantum dots are modeled as dipole emit- 
ters and the underlying problem is then that of a single dipole emitter in proximity to a 
MNP [H SI El El El [TOl HIl [12]. In general, when the influence of a particle on the dipole 
emission is studied, the outgoing electromagnetic fields of a dipole emitter are expanded 
into the electromagnetic multipole fields centered at a origin of and incident on the parti- 
cle [H El El [131 ttH US] ■ This transformation induces formally an infinite number of higher 
order multipole electromagnetic fields that can be characterized by the angular momen- 
tum numbers I and —l<m<l [U El El [131 [111 US]. A recent publication by Carminati et 
al [10] suggested an approximate analytical approach to calculate the fluorescence decay 
rate and the radiative and non-radiative rates. The rates were determined solely as a 
function of the particle polarizability a. For instance, the fluorescence decay rates for 
a parallel (tangential) and perpendicular (radial) dipole orientation with respect to the 
particle surface have been determined as (cf. Eqs. (9) and (10) of Ref. [TO] ) 



where r is the (radiative) decay rate in free space, z is the distance to the center of a 
nanoparticle and k stands for a wave vector. 



As in Ref. [10], let us focus on the distance dependence of the various rates of a dipole 

emitter at a proximity to Ag spherical nanoparticle (AgNP) with radius of 5 nm. Two 
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different emission wavelengths of the dipole emitter will be considered: 354 nm and 612 
nm. The dielectric constant of AgNP at the respective emission wavelengths is taken to 
be e(354 nm) = -2.03 + 0.6i and e(612 nm) = -15.04 + 1.02i pH [16]. For the sake of 
simplicity it will be assumed throughout that the host medium is air. Already a first glance 
at the fluorescence decay rates in Figs. [T] and [2] suffices to appreciate marked differences 
between the rates obtained using the approximate analytical approach of Carminati et 
al and the exact results [H El El [131 EH EH]- The converged exact rates in Figs. [1] and 
[2] were obtained using the freely available code in CHEWFS [17] (see Appendix [A] for a 
brief description of the code) and comprise contributions of the induced multipoles with 
I up to l max = 50. Note in passing that the rates labeled as "exact" in Figs. 1 and 2 of 
Ref. [10J are in fact those given by Eqs. (CQ). The latter correspond to a situation when 
all the induced higher order (I > 1) multipole contributions to the real exact decay rates 
but the dipole one (1 = 1) are neglected. This is confirmed by Figs. [1] and [21 wherein 
the nonconverged calculations of the exact decay rates with a cut-off at I = 1 are shown 
to coincide with the approximate rates. Although the authors of Ref. [10] warned at 
one occasion that the validity of their expressions (PQ) requires the distance z to remain 
larger than a few radii of the MNP, the warning remained largely ignored. Indeed, had 
the warning been taken seriously, any discussion of the apparent quantum yield as in 
Fig. 3 of Ref. [10] should have been avoided in the dipole approximation: to address the 
apparent quantum yield and the question of fluorescence quenching at short distances, 
one is required to calculate various rates for the very forbidden distances z within a few 
radii of the MNP. In this work we discuss in depth the following issues: 

1. The derivation of the particle polarizability a in [10J misses a dynamic depolarization 
term. 

2. Even if both the radiation reaction correction and dynamic depolarization are in- 
cluded in a, there is, for the examples studied in [10], essentially no difference with 
regard to the corresponding rates obtained by substituting the static Rayleigh po- 
larizability a.R in the formulas of [10] . 

3. For the examples studied in [10], the dipole-dipole coupling cannot describe satis- 
factorily the non-radiative rates even if the emitter is as far as three particle radii 
from the particle surface. 

4. The familiar Gersten and Nitzan quasistatic theory [18] yields non-radiative rates 
that agree much better with the exact rates than those of Ref. [10] . 

5. Because the validity of expressions ([1]) requires the distance z to remain larger than 
a few radii of the nanoparticle, the dipole approximation of Ref. [ID] cannot be 
applied to describe quenching [201 Ell 1221 [23] . The dipole-dipole approximation can 
only yield an upper bound on the apparent quantum yield. 

6. All the rates are affected by the size corrections to the metal dielectric function 

[21 [221 [25]. 
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The paper is organized as follows. Sec. [2] introduces a hierarchy of approximations for 
the particle polarizability and studies the effect of a particular choice of the polarizability 
on the resulting rates. Although the results of Ref. [10J have just been shown to correspond 
to the nonconverged calculations of the exact decay rates with a cut-off at / = 1 (see 
Figs. [Hand [2]), the behaviour of the radiative rates will be shown to be different from 
that of the non-radiative rates in that the / = 1 approximation to the radiative rates 
is essentially exact in the studied situation. The latter point is the subject of Sec. [31 
The remaining points such as a comparison with the Gersten and Nitzan theory [18], the 
distance dependence of non-radiative rates, some points regarding the calculation of the 
power P a b s absorbed inside the particle, and the influence of the size correction to the bulk 
dielectric function on the rates are discussed in Sec. HI We then conclude with Sec. El A 
brief description of the freely available codes, which were used to generate the results of 
this paper, is provided in Appendix [A] 



2 A hierarchy of approximations for the particle po- 
larizability and corresponding rates 

In Ref. [TO] the particle polarizability was given by 

acM = ^5 7T' ( 2 ) 

with an (denoted as ao in Ref. pH]) being the familiar static Rayleigh polarizability of a 
sphere, 

where a and e(u) are the radius and the dielectric constant of the particle, respectively. 
The fc 3 -dependent term is the familiar radiative reaction correction emphasized by earlier 
Ref. [27]. The radiative reaction correction applies to any oscillating dipole, be it an 
elementary molecular dipole or the dipole induced on a small (nano)particle, and results 
from the Abraham- Lorentz equation (see, e.g., Sec. 16.2 and Eqs. (16.8-9) of Ref. [29J). 
However, the exact polarizability for a sphere is obtained from Mie's solution as 

67T 

OLMie = -ljpT E l, (4) 

where is the electric dipole T-matrix element (the latter corresponds to the minus of 
the electric dipole Mie coefficients as given by Bohren and Huffman [30]). Upon introduc- 
ing the size parameter x = ka, one finds in the long-wavelength limit [31] 

3 e ~ 1 

OLMie ~ OlMie;as = 47TO 



e + 2 - (6e - 12) % _ ^ ( £ _ i) 

(5) 
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To this end note that each of the polarizabilities can be parametrized as 



Here A 2 = A 2 (e) (A 2 = for a R and ac) is in general a function of e, and hence a complex 
function of frequency, whereas A 3 (A 3 = for a^; A 3 = 1 for etc and ctMie;as) is a rea/ nu- 
merical constant. The special case A 2 = A3 = 1 of the above parametrization corresponds 
to the approximate Meier and Wokaun polarizability ctuw [2H1 E]. The polarizabilities 
a R , ac, cymw, and ctMie;as in this order approximate the particle polarizability ctuie with 
an increasing precision. For instance, when ac is compared against otuw or auie-^s, one 
immediately notices a missing fc 2 -term in the denominator of ac- The missing term is 
the so-called dynamic depolarization term [2EJ EJJ. For a given constant homogeneous 
external field, the Meier and Wokaun polarizability (Xmw presumes a constant internal 
field within the particle. In contrast to a mw, ®-Mie;as takes implicitly into account an 
order x 2 deviation from the constant internal field |31j . 

TABLE 1. The values of various polarizabilities a =(Re a, Im a) of AgNP for 

A = 354 nm and A = 612 nm. 





A = 354 nm 


A = 612 nm 


a R /(4:7ia 3 ) 
a c /{4:7ta 3 ) 
olmwI (47ra 3 ) 

ocMie/ (47ra 3 ) 


(1.249, 4.988) 
(1.244, 4.977) 
(1.055, 5.068) 
(1.119, 5.036) 
(1.093, 5.042) 


(1.229, 0.0179) 
(1.229, 0.0180) 
(1.233, 0.0181) 
(1.231, 0.0181) 
(1.231, 0.0181) 



As shown in Ref . [31] , the dynamic depolarization term can also be derived from the very 
same Green's function approach of Ref. [10] , the fact that has eluded the authors of Ref. 
[TO] . Without the dynamic depolarization term, the polarizability ac does not take into 
account any size dependence of the dipolar localized surface plasmon resonance (LSPR) 
position. The size dependence is determined, up to the order of x 2 , by the following 
equation for the real part of e = e' + is" (see Sec. 12.1.1 of Ref. [30J): 

(7) 

Although the dynamic depolarization term would cause for a > 40 nm an appreciable 
difference of the rates obtained with aMie-,as or auie compared to those with ac, with the 
radius of AgNP of mere 5 nm and x < 0.09, there is, in the present case, essentially no 
difference as to what polarizability is substituted into the approximate rate equations of 
Ref. [10J. The latter point is illustrated by Table 1 that lists the values of various polariz- 
abilities, which have been introduced so far. Therefore, the difference in the fluorescence 
decay rates obtained using the approximate analytical approach of Carminati et al [10] 
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with regard to the exact results [U El El ESI EH E5] (see Figs. El and cannot be the 
result of a wrong choice of polarizability. 

To this end note that the implications of the unitary bound on the physical S matrix, 

< £7£7+ < 1. (8) 

Given the parametrization £7=1 + 2zT of the S matrix, the unitarity bound implies for 
the physical T matrix and the polarizability 

- 7 - \T\ 2 < Re T < -|T| 2 < 0, |cf < Im a < ^ \a\ 2 + (9) 

4 D7T D7T Ik 6 

where the second bound follows from the first one on substituting ik 3 a/(Q7r) for T. Because 
Im olr = for a real e, the familiar static Rayleigh polarizability a# violates the above 
physical bounds. The violation is reflected by the fact that a tot ^ a sca + a a b s , where a to t, 
o'sca, Cabs are the total, scattering, and absorption cross-sections. Nevertheless, since a to t 
is in the present case dominated by a a b s , one could in the present case calculate the rates 
with the familiar static Rayleigh polarizability without ever noticing any difference 
with respect to the results obtained by means of ac (see Figs. [3]to|6]). 



3 Radiative and non-radiative rates 



The non-radiative rates are derived from the power P a b s absorbed inside the particle 
[H El EDI EE] • Carminati et al pU] derived the following approximate expressions for the 
rates (see Eqs. (15) and (16) of Ref. [TO]): 
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(10) 



The respective radiative rates were then derived in Ref. [TU] as a difference of the cor- 
responding total decay rates in the short distance limit ka < kz <C 1 (Eqs. (11), (12) in 
Ref. [10]) and the non-radiative rates (fTOl) . The resulting approximate expressions for the 
radiative rates were (Eqs. (17) and (18) of Ref. [TO]): 
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It can be verified that the above formulae of Ref. [10] hold irrespective which polarizability 
of Sec. |2] is used. A necessary and sufficient condition is to show that the power P a b s 
absorbed inside the particle remains to be given by the formula (see Eq. (13) of Ref. [10] ) 
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The latter formula can easily be verified to follow from the more general expression 

p - U£ ° „ IF I 2 ft*\ 

r abs qJs, u o,bs\ I - J exc\ i V 1J / 

where a a b s is the absorption cross-section of a MNP. Indeed, provided that one factorizes 

the electric dipole element of the T-matrix (in the SI units) as T E1 = i^a [cf Eq. PJ], 
the total (extinction), scattering, and absorption cross sections can be expressed as 

67T 

Otot ~ -- rj Re t ei = k Im a, (14) 

F |T B1 | =-H, ^ (15) 

a a6s « ^ (1-|1 + 2T S1 | 2 )=A; (ima-^^ 2 ). (16) 



Consequently, upon substituting (|T6|) into (|T3j) one recovers the formula ( !T2|) . irrespective 
of the choice of the polarizability of Sec. |2j 

A comparison of the approximate and exact radiative rates is shown in Figs. |3] and HJ 
The results show clearly that all rates coalesce to a single line, irrespective if obtained in 
the dipole-dipole approximation with ocr, ac, ctMie;as, or a Mie , or as the exact radiative 
rates with a cut-off imposed at either I — 1 or I — 50. Thus the induced multipoles of the 
order I > 2 have in the present case hardly any effect on the radiative rates, because the 
exact radiative rates are very precisely approximated by taking into account solely the 
dipole-dipole interaction. The latter means that the conclusions of Carminati et al [10] 
regarding the distance dependence of the radiative decay rate are correct. The distance 
dependence is chiefly dominated by a z~ 3 dependence, with a z~ 6 dependence being visible 
at plasmon resonance. On the other hand the results are hardly surprising and appear to 
have been anticipated by Gersten and Nitzan, who also approximate the exact radiative 
rates by taking into account solely the dipole-dipole interaction (see Eqs. (B.18') and 
(B.43') of Mathematical Appendices to Ref. [Ig]). 

An analogous comparison for the non-radiative rates is shown in Figs. OandEl Because 
the leading contribution to the non-radiative rates ([TO]) is proportional to Im a, and a 
reaches a maximum at the LSPR, the non-radiative decay is strongly enhanced when 
the emitter radiates at the plasmon-resonance frequency of the nanoparticle. Again all 
the approximate dipole-dipole interaction results, irrespective if obtained with (Xr, ac, 
&Mie;as, or OLMie-, together with the exact radiative rates obtained with a cut-off imposed 
at I = 1 coalesce to a single line. However, contrary to the radiative rates shown in Figs. 
|3] and HI the non-radiative rates cannot be approximated by taking into account solely 
the dipole-dipole interaction. Clearly, as demonstrated by the exact non-radiative rates 
with the cut-offs imposed at I — 2, I = 4, and I = 50, respectively, the contribution 
of induced multipoles of the order / > 2 has a pronounced effect on the non-radiative 
rates. The latter shows up in the distance dependence of the non-radiative rate. This 
point is further elaborated in Sec. 14.21 Because the exact T nr is always higher than the 
approximate one, whereas the exact and approximate T r are essentially equal, the dipole- 
dipole approximation of Ref. [10J can only provide an upper bound on the apparent 
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quantum yield Y that is denned as the ratio Y = T r /(T r + r nr ). The different distance 
behavior of the radiative and non-radiative rates implies that the apparent quantum yield 
always vanishes at short distance, resulting in a quenching of the emission. 

The results for the exact radiative and non-radiative rates shown in the figures were 
obtained using the freely available code CHEW [32] (see Appendix[S]for a brief description 
of the code). The exact rates comprise contributions of the induced multipoles with / up to 
Imax = 50. The exact radiative rates are derived from the total radiated power. The latter 
is calculated by the code CHEW by means of Eq. (108) of [6], where, in the notation of Ref. 
[6], the functions f^i are given as linear combinations of spherical functions (Eqs. (69), 
(72) of [6]) taken at the dipole position, with the coordinate origin located at the sphere 
center. [In the special case of a dipole outside a homogeneous sphere, the coefficients Q 7 ; 
(Eq. (72) of [B]) that enter the definition of f 7 i (Eq. (69) of [BJ) reduce to the Mie's 
expansion coefficients (cf. Eqs. (32), (33) of (6]).] The exact power P a b s absorbed inside 
the particle is calculated by the code CHEW by means of Eq. (119) taken in combination 
with Eqs. (116) and (120) of [BJ. In the special case of a dipole outside a homogeneous 
sphere one has, in the notation of Ref. [BJ, N — n — 1, and the coefficients a 7 / and in 
the integrands of the respective integrals Af in Eq. (116) of [6] reduce to 

a,i = VW(2) = 1/T+ 7Z (1), 6 7/ = 0. (17) 

It is worth to emphasize that the exact rates shown in the figures were all calculated 
independently: (i) the fluorescence decay rates have been determined in CHEWFS from 
the imaginary part of a relevant dyadics (cf. Eq. (1) of Ref. [10] Eq. (135) of Ref. 
[6]), the radiative rates have been determined in CHEW from the total radiated power 
of a dipole escaping to infinity (Eqs. (108) and (124) of Ref. [6]), and the non-radiative 
rates have been determined in CHEW from the total Ohmic loss of dipole radiation (see 
Eqs. (Ill), (115), (116), (133) of Ref. [BJ). One can verify that the sum of the radiative 
and non-radiative rates as obtained using CHEW coincides with the total decay rate as 
obtained using CHEWFS. 

4 Discussion 

4.1 A comparison with the Gersten and Nitzan rates 

In the special case of a dipole outside a homogeneous sphere, the functional form of the 
respective exact expressions for the total radiated power [Eq. (108) of [6] with / 7 / = 
ji(krd) + Kryih\ ; (&r<i), where k 7 ; are the Mie's coefficients (ai for 7 = M and 6/ for 7 = M; 
cf Eqs. (69), (32), and (34) of [6])] and for the power P a b s absorbed inside the particle 
[Eq. (119) of [6 J with c? 7 / = h^\krd)} look quite similar. However, whereas in Eq. (108) 
of [6] the Hankel functions are with increasing I multiplied by rapidly decreasing Mie's 
coefficients, in Eq. (119) of [BJ the Hankel functions are multiplied with increasing / by 
much slowly decreasing values of the integrals iff (Eqs. (116) of [6J) taken over the 
particle volume. The latter explains, which is also an intrinsic feature of the approximate 
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quasistatic Gersten and Nitzan theory [T5], that in the small particle limit (i) I = 1 is 
enough to describe the radiative rates (see Eqs. (B.18') and (B.43') of Mathematical 
Appendices to Ref. [H]), whereas, in principle, (ii) an infinite series in I is required to 
describe the non-radiative rates (see Eqs. (B.24') and (B.45') of Mathematical Appendices 
to Ref. [H]). This is also the reason why the non-radiative rates calculated by the Gersten 
and Nitzan theory agree much better with the exact electrodynamic non-radiative rates 
as confirmed by Figs. [7J and In arriving at the results, the same formulas as Eqs. (3) 
and (4) of Ref. [12] were used with the respective terms |1 — A| 2 , where A(_L) and A(||) 
are the Gersten and Nitzan image enhancement factors, being essentially equal to one. (In 
the present case, the ratios A/r r ini , with T r>int denoting the intrinsic vacuum radiative 
rate, are of the order of 10~ 10 .) That is analogous to Ref. [19J, where |1 — A| 2 f» 1 is 
implicit upon comparing Eqs. (5) and (6) with the exact formulas (B.18') and (B.43') of 
Mathematical Appendices to Ref. [18J. The above conclusions, i.e. that the molecular 
dipole excites all the multipoles of the sphere, but only the dipole radiates, whereas higher- 
order multipole terms gives rise to dissipation, have been confirmed within the quasistatic 
approximation also by other researcher [3U [351 EH] . 

The Gersten and Nitzan results are supported by experiment [12]. A mismatch be- 
tween the theoretical and the measured non-radiative rates reported in earlier work [H [7] 
may partly be attributed to nonlocal effects [33J and, partly, to the fact that a point 
dipole may be a too crude approximation for fluorophores but remains rather good for 
phosphors (the latter have a dipole strength that is reduced by orders of magnitude com- 
pared to that of fluorophores [12]). The nonlocal effects, in general, lead to significantly 
greater fluorescence rates and smaller non-radiative decay rates for the admolecules [33] • 
Eventually, Sec. 14.41 below demonstrates that the theoretical values of non-radiative rates 
for fluorescence near a LSPR can be further decreased toward the measured non-radiative 
rates [U [7J by acounting for size corrected dielectric function. 

Although the approximate quasistatic Gersten and Nitzan theory [18] has been shown 
to be superior to the dipole-dipole approximation of Ref. [10], the use of the Gersten and 
Nitzan theory [TB] is also not without limitations [12] . For instance, in contrast to the 
exact Chew expressions for the radiative rates [131 E] , the theory of Gersten and Nitzan 
[18] does not yield the correct asymptotic results for a flat surface in the limit a — > oo 
with a fixed value of the separation d of the dipole from the surface of a sphere. Moreover, 
as emphasized by Kim et al [19], in addition to the conditions a <C A and rf C A, the 
value of sphere radius a cannot be too large for a fixed d if the Gersten and Nitzan model 
[T8] . and hence any dipole-dipole coupling model, is to be applied for the radiative rates. 
Eventually, contrary to the exact electrodynamic results, the ratios r nr /Fo and r r /Fo, 
depend on the intrinsic vacuum radiative rate, because both image enhancement factors 
A(_L) and A(||) of the Gersten and Nitzan model depend on T rjint [18]. 

4.2 The distance dependence of non-radiative rates 

The approximate non-radiative rates f flO|) gave an impression that their short distance 
dependence is, like in the Forster energy transfer, characterized by a z~ 6 -like dependence. 
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However, it has been just shown (see Figs. [5] and [6]) that a dipole-dipole coupling is not 
appropriate to describe the non-radiative rates of a dipole emitter in a proximity of a 
MNP. Taking for example r nr .||, the exact r nr .|| for the emission wavelength of 612 nm 
shown in Fig. [5] is at the distances of 5 nm, 10 nm, and 15 nm from the AgNP surface 
by 101%, 33%, and 17% higher than the approximate dipole-dipole contribution of Ref. 
[TO] . Similarly, the exact r nr; || for the emission wavelength of 354 nm shown in Fig. [6] is 
at the distances of 5 nm, 10 nm, and 15 nm from the AgNP surface by 56%, 20%, and 
11% higher than the approximate dipole-dipole contribution of Ref. [10J. Because of the 
importance of the higher-order (I > 1) multipoles contribution, which decay as z~( 2l+4: \ it 
makes rather sense to consider the dependence of the non-radiative rates on the separation 
d of the dipole from the surface of a sphere [3U [3?]. Depending on the particle radius, 
the non-radiative rates can decay between « d~ 5 and d~ 3 , wherein a near d~ 3 decay is 
already observed for the particle radii > 10 nm (see Figs. 16 and 17 in Sec. 3.6 of Ref. 
[31] and Fig. 2.15 of Ref. [3?]). A slower than d -6 dependence of the non-radiative rates 
is valid even for particles as small as 1 nm in diameter [37]. Upon taking into account the 
results by Dulkeith [37], one can conclude that even for MNPs with the radius as small as 
0.5 nm and the emitter as far as the distance of 3a from the particle surface, the emitter 
and the MNP do not interact exclusively through dipole-dipole coupling. 

A full quantum-mechanical calculation for a spherical metal surface obtained within 
the time-dependent local density approximation (TDLDA) by Ekardt and Penzar [38] 
essentially confirms the above picture. Formally similar to the classical electrodynamic 
case [U [6], the quantum- mechanical non-radiative rates (see Eq. (5) in Ref. [38]) were 
given as an infinite series over I of the terms proportional to (Im ai)/z 2l+A , where ol\ is the 
/-pole polarizability of the sphere. The difference with the classical case is in that oli{oj) 
is the quantum-mechanical dynamical /-pole polarizability calculated within the TDLDA. 
Because each term in the series has a different distance dependence, the overall distance 
dependence of the non-radiative rates cannot be uniquely determined, simply because 
everything is mixed up [38] . Note in passing that the TDLDA results incorporate the 
so-called nonlocal dielectric response [331 EHJ EE] without any approximation. 

Unfortunately, a detailed comparison of the classical electrodynamic results [TJ [6] 
against the full quantum-mechanical calculation remains to be an open question. An 
obstacle in comparing quantum-mechanical TDLDA with the classical electrodynamics is 
that the former works with an electron density parameter [38], from which there is still 
some distance to calculate the dielectric function of say silver, which would then serve as 
the input parameter for the classical electrodynamic calculations. Thus at present it is 
not possible to answer conclusively at which distances from a MNP surface the classical 
results can no longer be applied. 

4.3 The power P a & s absorbed inside the particle 

It turns out that it is more reliable to calculate the power P a b s absorbed inside the particle 
by using the formula f TTB"]) . On the other hand, a direct integration approach of Carminati 
et al [JD] leads to a contradiction. Indeed, according to the definition of the polarizability, 
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the field inside the particle, Ej n , is related to the external field E ea;c acting on the particle 
as follows: 

_ 3a ^ 
E • = E 



A-na 3 (e - 1) 

On using the parametrization ([6]) of the polarizability, 
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When this expression is substituted into the formula (B.l) of Ref. [TO 
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With the help of the identity (B.3) of Ref. [TO], 
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and upon comparing with (cf. Eq. (B.5) of Ref. [TO] 
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whereas Im A 2 = for or, and «mvk- Upon comparing the two expressions (TT21 
and fl24|) for the absorbed power P abs one finds that, in the case of aMie-,as, they differ 
in the term proportional to Im A 2 . The result fl24l) would imply a modification of the 
non-radiative and radiative rate formulae of Ref. [TO]. However, the latter difference of 
the two expressions for P abs appears to be a consequence of that Eq. (|2"U|) is only an 
approximation to 



P, 



abs 



\V m \ 2 dV, 
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where the integral is over the volume of the particle. Eqs. (1121) and (12T)]) are equivalent. 
However, the hypothesis that, for a given constant homogeneous external field, the internal 
field inside the particle is a constant, and which leads to Eq. (120]). introduces an order x 2 
deviation from the exact result. Although at a first approximation the internal field inside 
the particle can indeed be assumed to be a constant, this is of course not true exactly 
|41j . It has been shown that the very form of the dynamic depolarization term in ctMie;as 
implicitly implies an order x 2 deviations from the constant internal field [31]. The order 
x 2 deviation from the exact result for P a b s when using Eq. ( 120]) is of the same order as 
the Im A 2 term that makes the difference between Eqs. ( fl~2]) and (124]) . It is therefore 
natural to assume that if the volume integration in Eq. ( 126]) were performed by taking 
into account the order x 2 deviation from the constant internal field, one would arrive at 
the result ( I12p obtained directly from the absorption cross section. 

4.4 Effect of size corrections to the bulk dielectric function 

It is generally accepted that for radii of MNPs smaller than the mean free path of con- 
duction electrons (42 nm for Au and 52 nm for Ag) the dielectric function of free electron 
metals should be corrected due to a reduced mean-free path of the electrons because of 
the scattering of the electrons on the MNP surface [23] 

e s (u) = e b (u)-e°(uj)+e°(uj). (27) 

Here, e b is the bulk metal dielectric function, and is the Drude dielectric function 

e^u) = l-u 2 p /[uj{uj + i 1 )l (28) 

describing the bulk metal conduction electrons with the plasma frequency u p and the 
bulk damping constant 7. The term ef d differs from in that 7 is replaced by the 
size-corrected value [23] E6] 

7sd = 7 + Av F /a, (29) 

where vf is the Fermi velocity and A (not to be confused with A'jS in the preceding sec- 
tions) is a phenomenological fitting parameter [24J. For gold and silver the size-correction 
becomes appreciable when, as in the present case, the sphere radius is < 8 nm (see for 
instance Fig. 4 [25]), i.e. less than one fifth of the mean free path of conduction elec- 
trons. The following parameters were employed for silver: vf = 1-39 nm fs -1 , u p = 72700 
cm -1 , and 7 = 145 cm -1 [32] [33]. Given the above parameters, e(354 nm) = —2.03 + 0.6i 
would change into e(354 nm) ps —2.028 + 0.686i for A = 0.25 and into e(354 nm) ps 
-2.008 + 0.945i for A = 1. Similarly, e(612 nm) = -15.04 + 1.02z would change into 
^(612 nm) = -15.022 + 1.466i for A = 0.25 and into e(612 nm) = -14.849 + 2.789i for 
A — 1. Following experimental results for individual gold NPs down to a = 2.5 nm by 
Berciaud et al [25], the value of A = 0.25 could be employed on the grounds of the simi- 
larity in the band structure of silver and gold [25] . Depending on the sample preparation, 
the value of A = 1 may be suitable in other cases |24j . 
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As shown in Figs. [7] and El by a dash-dot line, the effect of the size corrected dielectric 
function on the non-radiative rates is rather pronounced for A = 1. At the wavelength of 
354 nm, the size correction has the effect of decreasing the non-radiative rates, whereas 
for the wavelength of 612 nm the size correction has the opposite effect. Although the 
above results may at a first glance appear counterintuitive, they can be explained rather 
straightforwardly. Take for instance the leading I = 1 contribution, which is well ap- 
proximated by Eqs. f lTUj) . Eqs. f lTUj) show that the leading I = 1 contribution to the 
non-radiative decay is controlled through Im a rather than e" . According to Eq. fl22|) . 
Im an depends on e" as 3/(e") 2 near the dipolar LSPR. Therefore, at the proximity of the 
LSPR, increasing e" results in smaller values of Im an that causes the non-radiative decay 
to decrease. Thus by properly acounting for the size corrected dielectric function at the 
proximity of the LSPR may yield an additional mechanism for decreasing the theoretical 
values of non-radiative rates toward the measured non-radiative rates [H E]- 

On the other hand, e' is typically ten times larger in magnitude than e" far away 
from the LSPR. Consequently a change in e" has only a minor effect on the denominator 
Eq. (f22l and Im an becomes essentially proportional to 3e". Thus far away from the 
LSPR increasing e" results in larger value of Im an that causes the non-radiative decay 
to increase. 



5 Conclusions 

The methods based on dipole-dipole coupling PJJ [J_B] have known limitations in describing 
the radiative rates [J_H] • In line with earlier quasistatic results [TBI IMl E3 EE] , the present 
work reaffirms that the dipole-dipole coupling severely underestimate the non-radiative 
rates even in the parameter range, where it rather satisfactorily describes the radiative 
rates. An approximate analytical dipole-dipole coupling approach, such as that of Ref. 
|10j . may fail without any warning even for relatively small particles. The dipole-dipole 
interaction results of Carminati et al [TO] cannot be improved by any generalization of their 
expression for the dipolar polarizability, either by supplying the dynamic depolarization 
term or by taking the full Mie's expression. To correctly describe the radiative and 
non-radiative decay of fluorophores, chromophores, and quantum dots close to a metallic 
nanoparticle, the use of a complete theory is recommended [U El El [131 EEU EE]- Taking a 
Ag spherical nanoparticle (AgNP) with the radius of 5 nm as an example, a significant 
contribution (between 50% and 101% of the total value) of the higher-order (/ > 1) 
multipoles to the non-radiative rates was found at the emitter distance of 5 nm from 
the AgNP surface. Consequently, a dipole-dipole approximation can yield only an upper 
bound on the apparent quantum yield. In general, the higher-order (/ > 1) multipole 
contributions to the non-radiative rates cannot be neglected, even if a particle has a radius 
as small as a = 0.5 nm and the dipole emitter is further than 15 nm from the particle 
surface. The effect of the higher-order multipole contributions shows up in a markedly 
slower than d~ 6 distance dependence of the non-radiative rates. The non-radiative rates 
calculated with the approximate quasistatic Gersten and Nitzan method, which takes into 
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account the contribution of higher-order (/ > 1) multipoles, were found to be in much 
better agreement with exact electrodynamic results. Finally, the size corrected metal 
dielectric function was shown to decrease the non-radiative rates near the dipolar surface 
plasmon resonance. 
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A A brief description of the codes CHEW and CHEWFS 
with instructions how to adapt them to case of a 
homogeneous sphere 

The Fortran F77 codes CHEWFS [17] and CHEW [32] , which have been written according 
to the theory developed in Refs. [HI E] and which have been free available since October 
2004, are designed to calculate various decay rates of a dipole emitter interacting with a 
multi-coated sphere with an arbitrary number of shells. The emitter can be located either 
inside or outside a sphere. 

CHEWFS [IT] calculates the total decay rate and frequency shift directly from the 
respective imaginary and real parts of the Green's function at coinciding arguments (cf. 
Eqs. (135) and (137) of Ref. [6J). CHEWFS returns files dpfshiftm.dat (containing 
frequency shifts) and dipqmratesm.dat (containing total decay rates). 

CHEW [32] calculates the radiative and non-radiative rates, which are both calculated 
classically. The radiative rate is determined from the total radiated power of a dipole 
escaping to infinity. The latter is found by integrating the radial component of the time- 
averaged Poynting vector over the surface of a sphere with a radius going to infinity 
(cf. Eq. (88) of Ref. [6]). The output is in files dipratesm.dat and dipratesv.dat. The 
non-radiative rates are determined by CHEW from the Ohmic losses inside absorptive 
shells (cf. Eqs. (110) and (111) of Ref. [6]). The output is in files dipnrratesm.dat and 
dipnrratesv.dat, which contains non-radiative rates normalized to the radiative rates in a 
corresponding infinite space. 

The output files *m.dat and *v.dat differ only in that the rates are normalized with 
respect to those in the local medium at the dye position and in vacuum respectively. Each 
output file has data organized in four columns: normalized radial position (in the units 
of the bead radius) and the corresponding quantities for the respective tangential, radial, 
and averaged dipole orientations (from left to right). [The fourth column is not a simple 
average of the second and third columns but is calculated according to Eq. (123) of Ref. 
[6] as (r_i_ + 2T || ) /3)] . To modify the default versions of the two codes to the present case 
of a homogeneous sphere, the following steps have to be performed: 
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• The layers of a coated sphere with N shells are counted from one for the sphere core 
up to LCS=N+1 for the outermost shell and N+2 for the host medium. See also 
Fig. 1 of [6]. A homogeneous sphere has no shell and thus N=0. Consequently, 
the value of the integer parameter LCS, which controls the number of coatings, has 
to be set to LCS=1. Next the integer parameter ILCS has to be set to ILCS=1. 
The latter ensures that the dielectric constant specified by the complex parameter 
CCEPS becomes the dielectric constant of the sphere (core). 

• The value of the complex parameter CCEPS has to become one of the two values of 
the complex dielectric constant. At the same time the value of the integer parame- 
ter NMAT=0 has to be adopted. The latter choice corresponds to a dispersionless 
medium and should also be imposed if the data specified by CCEPS is given pref- 
erence. Otherwise, for NMAT > 1, the resulting value of the dielectric constant, 
which enters the Mie coefficients as the complex variable ZEPS1, is overwritten by 
either the Drude data or the real material data. 

• The integer parameter NABS, which specifies the total number of absorbing layers, 
and IABS, which carries the number of the absorbing layer, have to be both set to 
one. 

• The value of the real parameter NX specifies in the units of the sphere radius the 
length of the radial interval on which the decay rates are calculated. The value of 
NX is used to control the "do 200 il=l,nx" loop. 

• The angular momentum number / of the induced multipoles involved in the calcula- 
tion is controlled through the value of the variable LMX. The data setting LMX/1/ 
then forces the code to take into account only the dipolar contribution. The arrays 
of the codes are dimensioned with the integer parameter LMAX, where at present 
LMAX=60. The variable LMX can take any value between 1 and LMAX. 

B Supporting Information Available 

A MS Windows executable jH] is provided which calculates all the approximate rates 
shown here. 
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Figure captions 



Figure 1 - Normalized fluorescence decay rates versus the distance z to the center of a 
AgNP of radius a = 5 nm for the emission wavelength A = 612 nm. The exact rates shown 
by a solid black line here and in other figures of the paper comprise the contribution of 
the induced multipoles with I up to l max = 50. The exact rates with an imposed cut-off at 
/ = 1 are shown by a dash (red) line. The approximate dipole-dipole interaction results, 
irrespective if obtained with «c (dash-dot line; blue), aMie-, as (dot line line; green), or 
otMie (dash-dot-dot line; cyan), merge to a single line with that for the exact rates with a 
cut-off imposed at I = 1. 

Figure 2 - The same as in Fig. [T]but for the emission wavelength A = 354 nm. 

Figure 3 - Normalized radiative rates versus the distance z to the center of a AgNP of 
radius a = 5 nm for the emission wavelength A = 612 nm. The approximate dipole-dipole 
interaction results, irrespective if obtained with (short dash line; mangenta), ac (dash- 
dot line; blue), auie-^s (dot line line; green), or a Afe (dash-dot-dot line; cyan), merge to 
a single line with that for the exact radiative rates. There is essentially no difference if, 
in the latter case, the cut-off is imposed at either I = 1 (dash line; red) or I = 50 (solid 
line; black). 

Figure 4 - The same as in Fig. [3] but for the emission wavelength A = 354 nm. 

Figure 5 - Normalized non-radiative rates versus the distance z to the center of a AgNP 
of radius a = 5 nm for the emission wavelength A = 612 nm. The approximate dipole- 
dipole interaction results, irrespective if obtained with a^ (short dash line; mangenta), ac 
(dash-dot line; blue), auie-^s (dot line; green), or auie (dash-dot-dot line; cyan), merge 
to a single line with that for the exact radiative rates with a cut-off imposed at I — 1 
(dash line; red). The figure demonstrates a clear difference between the approximate 
non-radiative rates and exact non-radiative rates with a cut-off imposed at I = 2 (short 
dot line; purple), I = 4 (short dash-dot line; dark cyan), and I = 50 (solid line; black), 
respectively. 

Figure 6 - Normalized non-radiative rates versus the distance z to the center of a AgNP 
of radius a = 5 nm for the emission wavelength A = 354 nm. The exact rates with an 
imposed cut-off at I = 1 are shown by a dash (red) line. 

Figure 7 - Comparison of the normalized Gersten and Nitzan non-radiative rates against 
the exact rates for a AgNP of radius a = 5 nm for the emission wavelength A = 612 nm. 
The exact results with a size corrected dielectric function of silver for A = 1 are shown 
by dash-dot line (olive). All the rates were calculated with a cut-off l max = 50. 

Figure 8 - Comparison of the normalized Gersten and Nitzan non-radiative rates against 
the exact rates for a AgNP of radius a = 5 nm for the emission wavelength A = 354 nm. 
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The exact results with a size corrected dielectric function of silver for A = 1 are shown 
by dash-dot line (olive). All the rates were calculated with a cut-off l max = 50. 
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